Atomic scale friction between clean graphite surfaces 
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We investigate atomic scale friction between clean graphite surfaces by using molecular dynamics. The 
simulation reproduces atomic scale stick-slip motion and low frictional coefficient, both of which are observed 
in experiments using frictional force microscope. It is made clear that the microscopic origin of low frictional 
coefficients of graphite lies on the honeycomb structure in each layer, not only on the weak interlayer interaction 
as believed so far. 

PACS numbers: 07.79.Sp; 46.30.Pa; 31.15.Qg 

Keywords: Atomic scale friction, Graphite, Molecular dynamics simulation, Frictional force microscope 



Friction is one of the most familiar physical phenomena 
and has been investigated from long ago due to its impor- 
tance in various kinds of machinery and for the understand- 
ing of dynamics of many systems in science. Last two 
decades, atomic scale friction has been attracted much at- 
tention for the understanding of fundamental mechanisms of 
macroscopic friction and for many fields of high precision 
engineering such as nanomachine| l). Frictional force mi- 
croscope (FFM) has played an important role in study of 
atomic scale frictionjll 0|. Behavior of frictional force in 
FFM experiments has been studied theoretically and numer- 
ically based on the model which consists of a single atom tip 
and the potential by a substrate surface |3[ |j, |5|, |fj |8|] . The 
magnitude of load in FFM experiments in layered materi- 
als, however, is much higher than that estimated numerically 
and theoretically 01- Some groups pointed out that a flake 
cleaved from the substrate plays a crucial role in FFM exper- 
iments in such systems 1 2, EII9I [Tol [Till . It is not appropriate 
to discuss frictional phenomena with a flake by using the sin- 
gle atom tip model. We investigate in this letter atomic scale 
friction between clean graphite surfaces by molecular dynam- 
ics (MD). The model employed here simulates friction in the 
FFM experiments on a highly oriented pyrolytic graphite sub- 
strate with a flake. The present work enables us to understand 
results of the FFM experiments more deeply. Graphite is one 
of the most important solid lubricants. It turns out that the 
microscopic origin of lubrication property of graphite lies on 
the honeycomb structure in each layer, not only on the weak 
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FIG. 1: (a)The schematic illustration of the model and (b) the flake 
and the substrate in the present study. The solid and dashed lines in 
(b) denote bonds in the flake and those in the substrate, respectively. 

interlayer interaction. 

The model consists of a monolayer graphite substrate, a 
monolayer graphite flake and a spring which drives the flake 
as shown in fig.^a). Atoms in the flake obey the following 
eqs. of motion, 
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Here txiq is the mass of the carbon atom, xf, a = x, y, z, are 
components of the i-th atom position vector and a dot above 
xf denotes the time derivative. The x and y-directions are 
shown in fig.nb) and the z-direction is normal to the surface 
plane of the substrate. Atoms in the substrate are fixed be- 
cause their deformations are negligibly small in the present 
parameter range. The first term in the r. h. s. of eq. ([0 is a 
damping term proportional to the velocity, which reproduces 
the energy dissipation to the substrate. We set a damping con- 
stant 7 — 1.14 x 10~ 4 nNs/m. This value of 7 holds the 
stability of the flake for the time step of MD simulation At 
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= 2.74 xl0~ 15 sec. The second term denotes the force com- 
ing from interatomic interactions between the flake and the 
substrate. Vsub denotes the sum of pair atomic interaction po- 
tential between an atom in the flake and that in the substrate. 
We employ the Lennard- Jones (12,6) potential with a cut off 
length R, outside of which the interatomic potential is negligi- 
ble. We set Lennard-Jones parameters and the cut-off length 
ase = 0.965 x 1(T 2 eV, a = 0.340 nm and R = 2.0 nm. 
The above values of e and a are often employed in the inter- 
layer atomic potential in bulk graphite 1 12]. The third term 
denotes the force by the intraflake atomic interactions. We 
employ linear spring potentials for the bond length and the 
angle between bonds of nearest neighbor carbon atoms and 
z-displacement of an atom against nearest atoms. They are 
employed in the study of lattice vibrations and specific heat of 
graphite II 311 and FFM images |5j|6J]. The fourth term denotes 
the force by the driving spring with the components of spring 
constants, k a , a = x, y, z. The spring corresponds to a tip 
and a cantilever in FFM. Nf, x^, and xg denote the number of 
carbon atoms in the flake, a-components of a center position 
of the flake mass and a base position of the driving spring, 
respectively. Here we set k x = k y = 1.5 eV/nm 2 , k z = 5.0 
eV/nm 2 . Flake size is not known experimentally. The sizes of 
the flake N{ in the present calculation are 6, 20, 42 and 110. 
The results shown below are those for N { = 110| 14]. The 
shape of the flake is shown in fig. 0b). The size dependence 
is discussed later, x-q obeys the following eq. of motion, 

R c m c x^ = -7±b ~ L + k z (x z c - Xg) . 

±£ = v?, ^ = Kf, (2) 

Here Rq, the ratio of the effective mass of the tip and the 
cantilever to the mass of the carbon atom, is assumed to be 
100.0. L and (V^ , Vg) denote the loading force and the lateral 
driving velocity, respectively. The values of these parameters 
are kept constant during driving. 

Equations of motion, Q and (0 are solved numerically by 
using the Runge-Kutta method. In the initial state the system 
is stable with the applied loading force L and no lateral ex- 
pansions of the driving spring. The configuration corresponds 
to the AB stacking structure of bulk graphite as shown in 
fig.Qb). The frictional force is defined as a driving direction 
component of the force of the driving spring, k a (xg — a;g). 

It is to be noted that deformations of the flake are negligi- 
bly small in the present calculation. This is because that the 
interlayer atomic interactions between the flake and the sub- 
strate under load are much weaker than the intralayer atomic 
interactions in the flake. It is also to be noted that the rota- 
tion of the flake around Xq in the x-y plane is inhibited due to 
the potential barrier against it. Due to these two features the 
dynamics of the center of the flake mass, xq, governs the dy- 
namics of the system. At first we focus on the typical motions 
of sc. Figure 13 shows the driving direction component of the 
center of the flake mass, xq, as a function of that of the base 
position of the driving spring, x-q, for L — 100 nN, Vs = 0.22 
m/s in the case driven along the x (a) and y (b) directions, re- 
spectively. The stick-slip motions of Xq are observed. These 
periods agree with the periods of the graphite lattice in each 
direction, 0.426 and 0.246 nm, respectively. 
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FIG. 2: xc as a function of xb driven along the (a) x and (b) y- 
directions for L = 100 nN, Vs = 0.22 m/s. The inset in (a) magnifies 
the behavior of sticking position D. The sticking positions A, D and 
G correspond to those in fig. [3] respectively. 
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FIG. 3: The trajectories of xc for L = 100 nN, Vs = 0.22 m/s. The 
solid and dashed lines indicate the trajectories of xc driven along 
the x and y-directions, respectively. The heavy and thin dotted lines 
indicate the equipotential lines of the substrate potential as a function 
of a?c under the condition that the separation between the flake and 
the substrate is set to the constant value, 0.273 nm, which is equal 
to the time averaged value of the separation during driving for L = 
100 nN. The unit of numbers on heavy dotted equipotential lines is 
eV. The inset shows the configuration of the flake and the substrate 
which makes the AB stacking structure of bulk graphite. The solid 
and dashed lines in the inset denote the bonds in the flake and the 
substrate, respectively. The shaded region denotes potential valleys. 

Figure[3]shows trajectories of xq in the x-y plane for L = 
100 nN, Vs = 0.22 m/s. The contour lines denote equipoten- 
tial lines of the substrate potential. The motion of xq starts 
from the initial position A. The solid line in the figure denotes 
the trajectory of xc driven along the x-direction. By apply- 
ing driving force xc goes slowly to the position B, which is 
the point of inflection of the total potential. As soon as xq 
arrives at B, xq slips to the next sticking position close to D 
through the position C, both of which are the positions of the 
substrate potential minimum. The sticking time around D is 
much shorter than that around A due to large elongation of the 
driving spring along the ^-direction and much longer than the 
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slipping time. The sticking times around A and D are about 
99.88 and 0.08 % of the total period of the stick-slip motion, 
respectively. As soon as xq arrives at E, which is also the 
point of inflection, xq slips again close to I, which is equiva- 
lent to A. Then xq repeats the periodic stick-slip motion. The 
dashed line shown in the figure denotes the trajectory of xq 
in the case driven along the y-direction. By applying driving 
force xc goes slowly to the position F, which is the point of 
inflection of the total potential. As soon as a?c arrives at F, a?c 
slips to the next sticking position close to G, which is the posi- 
tion of the total potential minimum. The sticking time around 
G is shorter than that around A due to elongation of the driving 
spring along the x-direction. The sticking times around A and 
G are about 70.0 and 30.0 % of the period of the stick-slip mo- 
tion, respectively. As soon as a?c arrives at H, the flake slips 
again to the next sticking position close to J, which is equiva- 
lent to A. Then xq repeats the periodic stick-slip motion and 
the trajectory of xq becomes zigzag. The slipping times are 
of order of ten pico-seconds. The slips of xq occur between 
positions close to the substrate potential minimum, in which 
the flake makes the AB stacking structure of the graphite, as 
shown in the inset of fig.|3] The atomic scale stick-slip motion 
results from the binding of the flake close to the AB stacking 
positions of the graphite flake on the graphite substrate. 

Figure |4] shows the frictional force F as a function of xr 
for Vs = 0.22 m/s driven along the x (a) and y (b) directions, 
respectively. F initially increases linearly with x-q due to the 
elastic deformation of the driving spring. After F takes the 
maximum value the curves exhibit periodic saw-tooth shapes, 
which are observed in experiments I2ll4l. ll5l.llal . due to the 
periodic stick-slip motions of the flake. The regions A, D and 
G as shown in the insets in fig.|4](a,b) correspond to sticking 
states around those in fig.|3] respectively. 

We define the maximum static frictional force F s and the 
kinetic frictional force Fk as the maximum frictional force 
during driving and a time averaged value of the frictional force 
during one cycle of the stick-slip motion, respectively. Fig- 
ure |3 shows F s and Fk as a function of the loading force, L, 
driven along the x and {/-directions. The dependence of F s 
and Fk on Vs is weak. This is because that even in the present 
calculation the time scale of the driving is much longer than 
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that of atomic one. The values of the frictional force in fig.|3] 
are obtained by extrapolating data for Vs = 0.22, 0.44 and 
0.88 m/s to Vs = m/s for the comparison with experiments, 
in which the typical magnitude of Vs is about 40 nm/s. F s 
and Fk linearly depend on the loading force and have a finite 
adhesive term a, which is the frictional force at L = nN. 
The frictional coefficient fi is defined as, 



F = fiL + a. 



(3) 



Here F denotes F s or Fk. The magnitude of the frictional 
force in the range of the loading force in fig. [5] is very low 
with the frictional coefficients fi s = fik — 0.013 driven along 
the x-direction and fi s — fXk = 0.0055 driven along the y- 
directions. Here /i s and fXk denote the frictional coefficients 
of F s and Fk, respectively. These values are close to those 
observed in the experiment, fi = 0.012 |2]. 

Here we discuss the mechanism of the low frictional co- 
efficients of graphite. It has been believed that the low fric- 
tional coefficient of macroscopic graphite surfaces or graphite 
solid lubricants results from the weak interlayer interaction. 
We obtain, however, the frictional coefficient /i = 0.099 for 
a single carbon atom tip on the graphite substrate for the 
same pressure value with the case of the flake, which is much 
larger than those in our simulation of the flake and that in the 
experiment^]. The difference of the frictional coefficients in- 
dicates an another mechanism of the low frictional coefficients 
between graphite surfaces. In the flake there are two kinds of 
lattice sites A and B catching different forces from the sub- 
strate as shown in fig.|6] The substrate potential for the flake 
K u b is the sum of the substrate potential for a single atom at 
the sites A, V s ^ b , and that at the sites B, V® h . That is, V sub = 

^iKub + -WfiKub' where -^A and denote the number of 
the sites A and B. In the flake = iVg. Here we consider 
the case that a?c moves straight between the nearest substrate 
potential minimums A and C shown in fig. [3] The atoms at 
sites A move from the potential minimum to the maximum 
and those at sites B move from the maximum to the minimum 
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FIG. 4: The frictional force F as a function of ib for Vs = 0.22 m/s 
and L = 1, 10, 50 and 100 nN and in the case driven along the (a) x 
and (b) {/-directions. 



FIG. 5: Fk and F„, as a function of the loading force, x and y 
indicate the direction of the driving. 
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FIG. 6: The stacking structure of a flake on a substrate. The solid 
and dashed lines denote bonds in the flake and those in the substrate. 
The sites of carbon atoms A and B catch different forces of the in- 
teratomic interactions between the flake and the substrate. The con- 
figurations A and C correspond with the stacking structures at the 
positions A and C in fig. [3] 
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FIG. 7: The substrate potential during the motion of the flake A — > C 
shown in fig. [5] The solid, dashed and dotted lines denote the poten- 
tial value per atom for the flake as a function of a moving-direction 
displacement of xc, Axe, and the potential value as a function of 
Axe for a single atom at the sites A and B, respectively. 



as shown in fig.0 During the motion the increase of V s ^ h and 
the decrease of V® h cancel out each other. Then the substrate 



potential for the flake V au b does not vary in the route between 
the nearest potential minimums. As a result the flat potential 
valleys appears for the flake in the shaded region shown in 
fig- El The valleys enable the flake to move easily and yield 
low frictional force. Due to this mechanism graphite becomes 
a good lubricant with low frictional coefficient. On the other 
hand the single carbon atom above the graphite substrate has 
large frictional coefficient because of the absence of such can- 
cellation mechanism of the substrate potential. 

The frictional coefficient mainly depends on the trajectory 
of xc- When the flake is large enough that the interaction be- 
tween the flake and the substrate is so strong as to overcome 
the stiffness of the driving spring, the dominant contribution 
to the acceleration of the flake is flake size independent value, 

(^A^mb + ^BKub)/2' 7l C> f° r the same separation between 
the flake and the substrate. Here 



denote the coordinate 



derivatives at the sites A and B, respectively. Then because 
the trajectory is independent of the flake size, the frictional 
coefficient is also independent of the flake size. The flake in 
this simulation is sufficiently large, so that the frictional co- 
efficient is independent of the flake size. In fact the differ- 
ence between the frictional coefficients for Nf = 42 and 110 
is much small as about 2 %. 

In the present work we have investigated atomic scale fric- 
tion between clean graphite surfaces by numerical simulation. 
The simulation reproduces atomic scale stick-slip motion and 
low frictional coefficients of graphite. The former is due to 
the binding of the flake close to the stacking configurations on 
the substrate. The latter results from the cancellation of the 
forces coming from the substrate between two kinds of lattice 
sites in the flake, which is discovered in this letter for the first 
time within our knowledge. 
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